

 ***graphs
 set scheme plotplain 

 
 
 /// Full sample overview (2.3.1)
			use "Temp/Master", clear 
		  
			drop if year<1900
		   
		 
		  
			gen tax_non_trade_real_median = tax_non_trade_real 
			collapse (mean) tax_non_trade_real  ordinary_real (median) tax_non_trade_real_median , by(year)	
			line tax_non_trade_real tax_non_trade_real_median ordinary_real year	,  plotregion(margin(0)) bgcolor(white)  bgcolor(white)  ytitle("Real revenue (labor days per capita)") xtitle("Year" ) yla(,nogrid) xla(,nogrid)  ylabel(0(20)120)  xscale(r(1900,2020)) xsize(8)       xlabel(1900(10)2020)   legend(pos(11) ring(0) label(1 "Taxes  excl. trade & resources (mean)") label(2 "Taxes  excl. trade & resources (median)") label(3 "Ordinary revenues (mean)") ) 
		 
			graph export   "Output/Graphs/Appendix/APP_FIG_2_3_1.pdf", replace


			
			 
 /// Direct taxes only   (2.3.2)
			use "Temp/Master", clear 
		  
			drop if year<1900
		   
		 
		  
 			collapse (mean) direct_real  tax_non_trade_real   , by(year)	
			line direct_real  tax_non_trade_real	year ,  plotregion(margin(0)) bgcolor(white)  bgcolor(white)  ytitle("Real revenue (labor days per capita)") xtitle("Year" ) yla(,nogrid) xla(,nogrid)  ylabel(0(10)60)  xscale(r(1900,2020))    xsize(8)    xlabel(1900(10)2020)   legend(pos(11) ring(0) label(1 "Mean of direct taxes") label(2 "Mean of all taxes (excl. trade & resources)")  ) 
		 
			graph export   "Output/Graphs/Appendix/APP_FIG_2_3_2.pdf", replace
		  
   
  
 
/// Balanced sample forced labour  (2.3.3a)
			use "Temp/Master.dta" , clear  
			 
			 
			drop if year<1900

			collapse (mean)  ordinary_real total_forced_low total_forced_high  if expansion == 0, by(year) 
			line  total_forced_low total_forced_high ordinary_real  year, plotregion(margin(0)) bgcolor(white)  bgcolor(white)   ytitle( "Real revenue (labor days per capita)",axis(1)  )   yla(,nogrid) xla(,nogrid)  ylabel(0(20)120)     xscale(r(1900,2020))       xlabel(1900(10)2020)   legend(pos(11) ring(0) label(1 " Ordinary revenues  incl. labor tax (lower bound)") label(2 " Ordinary revenues  incl. labor tax (upper bound)") label(3 " Ordinary revenues (monetary taxes only)") ) 
			 
			graph export   "Output/Graphs/Appendix/APP_FIG_2_3_3a.pdf", replace

			
/// Full sample forced labour	(2.3.3b)


			use "Temp/Master.dta" , clear  
			 
			 
			drop if year<1900

			collapse (mean)  ordinary_real total_forced_low total_forced_high  , by(year) 
			line  total_forced_low total_forced_high ordinary_real  year, plotregion(margin(0)) bgcolor(white)  bgcolor(white)   ytitle( "Real revenue (labor days per capita)",axis(1)  )   yla(,nogrid) xla(,nogrid)  ylabel(0(20)120)    xscale(r(1900,2020))       xlabel(1900(10)2020)   legend(pos(11) ring(0) label(1 " Ordinary revenues  incl. labor tax (lower bound)") label(2 " Ordinary revenues  incl. labor tax (upper bound)") label(3 " Ordinary revenues (monetary taxes only)")  ) 
			 
			graph export   "Output/Graphs/Appendix/APP_FIG_2_3_3b.pdf", replace

			
			
